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ABSTRACT 

We present the results from a spectroscopic survey for intervening Mg II 
absorption in the spectra of 381 background QSOs conducted at the Multi- 
ple Mirror Telescope. This survey complements our earlier SDSS EDR Mg II 
survey, extending our results to lower redshift (z — 0.15) and weaker Mg II 
A2796 rest equivalent width (Wq 2796 ~ O.lA). We confirm two major results 
from that survey: the transition in the Wq 2796 distribution at Wq 2796 ~ 0.3A, 
and the H^ 2796 -dependent evolution of the incidence of systems. The nature of 
d 2 N/ dzdWQ 2796 is consistent with the idea that multiple physically-distinct com- 
ponents/processes contribute to the incidence of Mg II absorption systems in a 
WVdependent manner and evolve at different rates. A significant decrease in the 
total proper absorption cross section is detected in our MMT data for systems 
as weak as 1.0 A < Wq 279& < 1.5 A at z < 0.4. We discuss this W -dependent 
evolution in the context of the evolution of galaxy structures, processes including 
superwinds and interactions, and damped-Lya absorbers. We also consider the 
possibility that the observed redshift and W^ 2796 dependence of the incidence 
of absorption in spectroscopic surveys for low-ion/neutral gas results from the 
effects of dust-induced extinction. 

Subject headings: galaxies: evolution — galaxies: ISM — quasars: absorption 
lines 



1. Introduction 

Over the past 25 years, spectroscopic QSO surveys have established the statistics of 
intervening metal-line absorption systems, most notably those selected via C IV and Mg II 
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resonance transitions (e.g., Weyman et al. 1979; Lanzetta, Turnshek, & Wolfe 1987; Tytler 
et al. 1987; Sargent, Boksenberg & Steidel 1988; Steidel & Sargent 1992, hereafter SS92; 
Churchill, Rigby, Charlton, & Vogt 1999, hereafter CRCV99; Nestor, Turnshek, & Rao 2005a, 
hereafter NTR05). These systems identify structures based on their gas absorption cross 
section, and are widely believed to select galactic structures and extended gaseous galactic 
halos (e.g., Bergeron & Boisse, 1991; Steidel et. al, 2002; Churchill, Kacprzak & Steidel 
2005). The results of the early Mg II surveys, which we summarize in NTR05, indicated 
that the distribution of rest equivalent widths of the A2796 line (H 7 ^ 2796 ) can be described 
by a power law or an exponential, and that the incidence of the strongest systems seems to 
be decreasing with decreasing redshift. In NTR05 we presented a large 1,300 systems) 
Mg II survey using the Sloan Digital Sky Survey (SDSS) Early Data Release (EDR) QSO 
spectra. We demonstrated that the Wq 2796 distribution is fit very well by an exponential 
over the range 0.3A < W£ 2796 < 5A and that previous power law fits over predict the 
incidence of strong lines. However, extrapolating the NTR05 exponential to weak (0.02A 
< Wq 2796 < 0.3A) systems under predicts their incidence when compared to the results of 
CRCV99. This apparent transition is consistent with the idea that weak Mg II absorbers 
are in part comprised of a population that is physically distinct from intermediate/strong 
Mg II systems (e.g., Rigby, Charlton & Churchill 2002). Since the transition occurs at the 
boundary of two different surveys, however, it is important to investigate the possibility that 
this is simply an artifact of some unknown systematic difference between the surveys. 

NTR05 also investigated the incidence of Mg II systems as a function of redshift for 
different ranges of W^ 2796 . The results were consistent with the total proper cross section for 
absorption decreasing with decreasing redshift, especially at z < 1, such that the stronger 
lines evolve faster. The apparent evolution can be described by a steepening of the W^ 2796 
distribution with cosmic time. However, the data only reach z > 0.37 and the evolutionary 
signal has high significance only for very strong (Wq 2796 > 2.0A) systems. Thus, extending 
the measurement of dN/dz to lower redshift is a major goal of this work. Furthermore, the 
evolution of Mg II systems over cosmic time has a more general significance, since Mg II 
absorbers trace neutral hydrogen. Extending the statistics of Mg II absorption to lower 
redshift is therefore important for the study of the evolution of the H I content of the 
Universe as well. In fact, results from this survey have already been used to study the 
incidence of damped Lya systems at low redshift (Rao, Turnshek & Nestor, 2005). 

Accurate measurements of the incidence of Mg II absorption systems as a function of 
W$ 2796 and redshift for weak to "ultra-strong" (4A < Wq 2796 < 6k) systems and from 
high to low redshift, together with kinematic and relative abundance information from high- 
resolution spectroscopy (see, e.g., Ding, Charlton, & Churchill, 2005, and references therein) 
and composite-spectrum analysis (e.g., Nestor, Rao, Turnshek & Vanden Berk, 2003), are 
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allowing the ability to understand the physical nature and evolution of the structures selected 
by such systems. Here we present the results of a spectroscopic QSO survey conducted at the 
6.5m Multiple Mirror Telescope on Mount Hopkins, AZ. The goals of the survey are: 1) to 
provide a survey large enough to detect a significant number of intermediate/strong systems 
(0.3A< W^ 2796 < 2k) and sensitive enough to measure weak (0.1 A < W^ 2796 < 0.3k) 
systems in a single survey; and 2) extend the measurement of their incidence to lower redshift 
(z ?s 0.15). In §2 we describe the observations and our data reduction procedure. We present 
the results in §3. In §4, we discuss the relevance of our results to the problem of understanding 
the nature of low- ion metal-line absorbers, and present our conclusions in §5. 



2. Observations 

Observations took place over a 42 month period from January of 2001 to May of 2004 
at the 6.5m MMT on the summit of Mount Hopkins, AZ. Useful data were collected on 
17 nights spread over seven observing runs. The MMT spectrograph was operated using 
the blue-channel optical layout with the 800 grooves/mm grating, which corresponds to 
0.75A/pixel and a resolution of approximately 2.2A. Exposure times varied with the QSO 
magnitude and varying observing conditions, but were typically ten to fifteen minutes per ob- 
ject. Approximately 900 spectra of nearly 400 QSOs were collected. Quartz lamp exposures 
were used to correct for pixel to pixel sensitivity variances, and comparison lamp exposures 
(usually He-Ne-Ar) were used to wavelength calibrate the spectra. The wavelength coverage 
varied slightly for each observing run, but typically spanned ~ 3190A to either 5170A or 
ss 5465A. 

The resulting data set comprised a total of 381 useful QSO sightlines. Continua were 
fitted to the reduced QSO spectra and Mg II doublet candidates found, inspected, and 
measured in a manner similar to that described in NTR05. The same software were used, 
though with slight modifications to account for the differences in resolution and wavelength 
coverage between the SDSS and MMT data. 



3. Results 

The final MMT-survey Mg II absorber sample consists of a total of 140 doublets. They 
cover the redshift range of 0.1973 < z < 0.9265 and have W^ 279e values ranging from 0.128A 
to 3.165A. 

Figure 1 shows the distribution of W^ 2796 . The redshift-path coverage as a function of 
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Wq 2796 is also shown. The resulting dN/dW^ 2796 is shown as the histogram in Figure 2. 
The solid line is a maximum likelihood fit to the data having Wq 2796 > 0.3A of the form 
dN/dW* 2796 = $^e~^ (see NTR05) with W* = 0.509±0.047 and N* = 1.089±0.121. The 
red-dashed line is the fit to the low-z data from the SDSS EDR survey, with a mean redshift 
(zabs) = 0.655. The systems found in the MMT survey have (z a i, s ) = 0.589. The points 
represent data from CRCV99, with (z abs ) = 0.9. The distribution determined from the MMT 
data is consistent with both the SDSS EDR and the CRCV99 results. Although there are only 
26 systems in the two lowest W$ 2796 bins of the MMT survey ON/dW^ 2796 for W$ 2796 < 0.3A 
is in very good agreement with the CRCV99 results and significantly (1.4(7 and 2.4<r) above 
the extrapolation of the single-exponent fit to dN/dW^ 2796 for W<$ 2796 > 0.3A. Though the 
CRCV99 data has a higher mean redshift, weak Mg II systems exhibit the least evolution. 
Correcting the Wq 2796 < 0.3A points for evolution has only a minimal effect on Figure 2 
and makes the consistency more robust. Thus, the results from the MMT survey confirm 
the upturn in dN/dW^ 2796 below 0.3A that was first identified in NTR05. Figure 3 shows 
^mmt determined for systems with Wq 2796 > 0.5A, along with the values from the SDSS 
EDR survey. The W^ 2796 distribution determined from the MMT survey is consistent within 
the errors with the SDSS EDR survey results, considering the observed redshift evolution. 

The distribution of redshifts for Mg II absorption systems found in the survey is shown 
in Figure 4. Also shown is the sightline coverage for W™ in = l.oA, 0.6A, and 0.3A. The 
resulting dN/dz values as a function of look-back time in a (il\,£l M ,h) = (0.7,0.3,0.7) 
cosmology are shown in Figure 5. The lowest redshift bin corresponds to the range of 
look-back time not covered by the EDR survey, except for the panels with Wq 2796 > 2.0A, 
which did not have enough absorbers to allow for more than a single bin. Also shown are 
the respective data from the EDR and CRCV99 surveys. The dashed lines are the no- 
evolution curves normalized to the binned EDR/CRCV99 data, excluding the MMT data. 
The 0.367 < z < 0.956 MMT results are consistent with the other studies. Below z ~ 0.4, 
the MMT results for systems with W^ 2796 < 0.6 A and W^ 2796 > 1.5 A are consistent within 
the large uncertainties with no evolution (however, the EDR results give significant evidence 
that the incidence of the strongest systems does indeed decrease relative to the no-evolution 
prediction with decreasing redshift). Also, the MMT result for 1.0 A < Wq 2796 < 1.5 A at 
0.195 < z < 0.367 is significantly 1.5er) below the EDR-normalized no-evolution curve. 

NTR05 demonstrated a significant detection of evolution in the incidence of W^ 2796 > 2A 
systems at redshifts z < 1. Although our MMT sample is not large enough to confirm this 
result for large Wq 2796 , we do detect evolution at lower redshift (0.195 < z < 0.367) for 
systems with 1.0 A < W^ 2796 < 1.5 A. Furthermore, the values of W* derived from the MMT 
data (Figure 3) are consistent with the NTR05 results whereby the evolution of the incidence 
of Mg II systems is described by a steepening of dN/dW^ 2796 with decreasing redshift. While 
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W* was measured from the SDSS EDR data down to (z) ~ 0.7, the MMT data extend this 
result down to (z) ~ 0.4. 

4. Discussion 

The incidence d 2 N/dW dz of absorption systems is proportional to the product of the 
number density of systems at a given W and redshift times their average absorption cross 
section. For the Mg II absorbers considered in this work, W A2796 is more a measure of 
kinematics than of column density, since many of the kinematically-distinct A2796 lines 
comprising an absorber (which are not individually resolved in the EDR or MMT data) are 
typically saturated. This is particularly the case for intermediate/strong systems. We have 
shown that the total proper cross section of Mg II absorbers decreases with cosmic time at 
a W^ 2796 -dependent rate. While the reason for the differential evolution is not yet clear, 
multiple physically-distinct populations/processes contributing to the incidence of lines non- 
uniformly in Wq 2796 would provide a explanation if those populations/processes evolve at 
different rates. 

In NTR05 we discuss evidences for Mg II absorbers tracing a variety of physical systems. 
However, the possibility that the observed evolution is to some extent an artifact of an 
evolving bias due to dust-induced extinction must also be considered. We further discuss 
these scenarios below. 

4.1. Population-dependent Evolution 

According to Rigby, Charlton & Churchill (2002), the number of individual kinematic 
components ( "clouds" ) in a given Mg II absorber obeys a Poissonian distribution, except for 
a large excess of single-cloud systems which account for 2/3 of weak systems and have their 
origins in a population of objects and/or processes distinct from stronger absorbers. Indeed, 
an extrapolation of our exponential fit to dN/dW^ 2796 (§3) to systems with W A2796 < 0.3A 
accounts for « 30% of the systems with 0.02A < W A2796 < 0.3A predicted by the CRCV99 
results. Thus, it appears that a fraction of the "weak" absorbers are physically similar to 
intermediate-strength systems, while the majority (the single-cloud systems) have a different 
physical nature. While intermediate-strength systems show no evolution to within fairly 
strong limits over the large range 0.5 < z < 2 (Figure 5, also see NTR05 for discussion), 
there is evidence that the incidence of weak systems exhibits a significant decrease above 
z 1 (Lynch, Charlton & Kim, 2005). It should also be noted, however, that most single- 
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cloud systems have Wq 2796 < O.lA (the limit of this study). Thus, multi-cloud weak systems 
may be in part responsible for the upturn between O.lA < Wq 2796 < 0.3A (see also Masiero 
et al., 2005). These systems have been linked in part to dwarf galaxies (e.g., Zonak et al., 
2004). 

Imaging studies of intermediate/strong Mg II absorbers (Rao et al. 2003, Belfort- 
Mihalyi et. al, in prep) have shown these systems to be associated with galaxies that are 
drawn from a variety of morphological types. Also, the results of Kacprzak, Churchill & 
Steidel (2005) indicate that, while galaxy orientation shows no correlation with absorption 
properties, interactions are an important factor in determining the kinematic spread of the 
absorption in intermediate strength systems. Unfortunately, only a handful of systems have 
been studied in depth with complementary imaging/spectroscopy of the absorbing galaxy 
together with high resolution spectroscopy to reveal the absorption kinematics. Steidel at 
al. (2002) have shown that extended rotating disks contribute partially to the incidence 
in a subset of systems chosen to be associated with inclined spirals. However, the edge-on 
spiral studied by Ellison, Mallen-Ornelas, & Sawicki (2003) is inconsistent with extended 
disk absorption. They suggest superbubbles as the cause of the absorption, as postulated by 
Bond et al. (2001a). The few Wq 2796 > 2A systems for which high-resolution spectroscopy 
exists show kinematics consistent with the superwind picture (Bond et al., 2001b), though 
it has also been suggested that galaxy pairs/interactions may contribute to the incidence of 
such systems (see, e.g., Churchill et al., 2000). 

Though the strongest systems have been poorly studied due to their relatively small 
incidence, there is evidence indicating that the populations and processes associated with 
these systems differ from those associated with intermediate strength absorbers. Turnshek 
et al. (2005, see also Nestor et al. 2003) demonstrate that the average metallicity in Mg II 
absorbers is strongly correlated with absorption kinematics, which suggests that the strongest 
systems are associated with more massive and/or evolved galaxies. Furthermore, Menard et 
al. (in prep) find that the average (over a large number of absorbers) degree to which QSOs 
are reddened due to foreground Mg II absorbers is strongly W / ( ^ 2796 -dependent for strong 
systems. Scenarios involving both superwinds and interactions give a natural explanation 
for the increased metallicity and reddening with increasingly stronger systems. Interactions 
can both strip gas, causing large line-of-sight velocity spreads (i.e., large Wq 2796 values) and 
induce star formation which in turn enriches the gas. In the superwind picture, relatively 
highly enriched gas expelled in superwinds dominates in the strongest systems, whereas 
much of the gas responsible for most intermediate-strength absorbers has remained at large 
galactocentric distances (impact parameters) since it first condensed. 

As mentioned by Heckman (2002), the appearance of superwinds is intimately related 
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to the star formation rate (SFR) per unit area in a galaxy. While Lyman break galaxies 
(z ~ 3) and moderate-redshift (1.4 < z < 2.5) star-forming galaxies do show evidence of 
superwinds (Pettini et al., 2001; Shapley et al., 2003; Steidel et al., 2004), the SFR densities 
of ordinary local spirals (Kennicutt, 1998) are insufficient. The incidence of superwinds 
should therefore decrease with decreasing redshift, especially at z < 1 since the global SFR 
decreases below this value (Hopkins, 2004, and references therein). The evolution of the 
galaxy pair-fraction and merger rates is uncertain, however. Some studies have found little 
evolution over < z < 1 (Lin et al., 2004; Carlberg et al., 2000), while others have found 
strong evolution (Le Fevre et al. 2000). 

We have also begun an imaging study of individual "ultra-strong" systems, having 
Wq 2796 > 4A. Though the numbers are still quite small, early results suggest that ultra- 
strong Mg II absorbers may select galaxies that are preferentially bright compared to the field 
population or those selected by intermediate/strong Mg II absorbers (see Nestor, Turnshek 
& Rao, 2005b). 

It is unclear, however, how DLA systems fit into the overall picture. According to 
Rao, Turnshek & Nestor (2005) the likelihood of a Mg II absorber having a neutral hy- 
drogen column density above the DLA threshold (Nhi > 2 x 10 20 ) rises from 15% for 
Wq 2796 ~ 0.6A to w 65% for Wq 2796 ~ 3.0A. However, DLAs preferentially select neither 
bright galaxies, as the "ultra-strong" systems may, nor high-SFR systems (see Hopkins, Rao 
& Turnshek, 2005), which are required for superwinds. If the incidence of the DLA, su- 
perwind, and interacting/galaxy-pair populations are evolving at a different rate than that 
for the population/processes traced (at least in part) by the "classic" Mg II absorbers, this 
should contribute to the differential evolution seen in Mg II dN / dz. As the role of the DLA 
sub-population provides important clues for the study of low- ion/neutral gas, sorting this 
out should be an important goal of future work. 



4.2. Dust 

The differential evolution in dN/dz could also be affected by an evolving dust bias. 
Vladilo et al. (2005) claim that extinction from dust increases linearly with Zn column 
density, from small (A v < 0.02 mag) levels at N Zn < 10 12 to significant levels at higher 
column densities (A v m 0.15 mag at N Zn ps 10 12 8 ). The results of Turnshek et al. (2005) 
indicate that the average Zn column density increases to approximately this level for strong 
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(Wq 2796 ps 3A) Mg II systems. 4 Indeed, the Menard et al (in prep) results for E(B-V) versus 
Wq 2796 are roughly consistent with that predicted by Vladilo et al., considering the Turnshek 
et al. WQ 2796 -Nzn relation. Thus, the largest Nz n systems (which preferentially correspond 
to the largest Wq 2796 values) maybe be missed due to the background QSOs dropping out 
of magnitude limited and/or color selected surveys such as the SDSS. If this is indeed the 
case, we expect that the W^ 279% -Nzn relation will level off at large Wq 2796 values as the 
largest individual N Zn systems drop out of the survey. Also, we would expect a turnover 
in dN / 8Wq 279& at large Wq 2796 as these systems are preferentially missed. 5 Both of these 
tests will be possible with larger Mg II absorber catalogs formed from the full SDSS survey. 
We also note that although recent studies such as those emerging from the CORALS survey 
(e.g., Ellison et al., 2004) claim to find no dust biases, the size of their survey is not large 
enough to study the rare, strongest systems. 



5. Conclusions 

We have detected and measured 140 Mg II absorption systems in the spectra of 381 
QSOs obtained at the MMT Observatory. These data include systems with A2796 rest 
equivalent width in the range 0.1 A < W^ 2796 < 3.2 and extend our determination of the 
incidence of Mg II absorbers down to redshift z = 0.15. 

Coverage of Wq 2796 across the value Wq 2796 ~ 0.3A is important since the distribution 
dN / 8Wq 279& is represented very well by a single exponential function above Wq 2796 ~ 0.3 
A, but diverges from the fit with an excess of absorbers below this value. This effect was 
first noted in NTR05 by comparing our SDSS EDR survey data for W* 279e > 0.3A with 
data from CRCV99 having W^ 2796 < 0.3A. Here, we confirm the nature of the ON/dW^ 2796 
distribution in a single survey thereby removing the concern that the putative upturn was 
an artifact of unknown systematic differences in the two independent surveys. 

By extending the redshift coverage of our Mg II survey to lower redshift, we sample 
systems down to a look-back time of only ps 2 Gyrs, compared to ps 4 Gyrs for SDSS samples, 
and thereby increase the total look-back time covered by ~ 1/3. The low-redshift regime 
is of particular importance, especially for intermediate/strong systems, since the apparent 
evolution in their incidence is easily manifested only at low redshift. While NTR05 report 



4 Note, that since these are averages, some systems are expected to have larger Zn columns, even at smaller 
Wo A2796 values. 

5 However, the presence of a turnover, if found, could also be explained by a physical upper-limit to the 
velocity spread of Mg II absorbers. 
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evolution in systems with Wq 2796 > 2.0A at redshifts z < 1, we now detect evidence for 
evolution in systems with W^ 2796 > l.oA at redshifts z < 0.5. The evolution apparent in our 
MMT sample is consistent with the nature of the evolution described in NTR05, whereby 
the incidence decreases from the no-evolution prediction (i.e., the total proper absorption 
cross section decreases) at a rate dependent on Wq 2796 . 

It seems likely that multiple populations and processes contribute to the total cross 
section of Mg II absorption in a W / A2796 -dependent manner, and that the different evolu- 
tionary rates of these populations and processes lead to the W / ( ^ 2796 -dependent evolution in 
dN/dz. It is also possible, however, that biases due to dust-induced extinction give rise to 
the detected evolution. It will be necessary to determine the relative contribution of the var- 
ious physically-distinct components/processes to the incidence of Mg II systems, as well as 
obtain stronger constraints on the evolution of the incidence, in order to better understand 
low-ion/neutral gas absorption systems. Determining the relative contribution as a function 
of Wq 2796 will help us further understand the evolution in time of each type (and, perhaps, 
vice- versa). The results of our ongoing studies, such as the imaging of "ultra-strong" sys- 
tems and the analysis of the full SDSS sample, are two such projects that should add to our 
understanding of the nature and evolution of these systems. 

We acknowledge support for this work from the NSF. We also acknowledge Brice Menard 
for helpful discussions and ongoing collaborations related to this work, and thank the MMT 
staff and the members of the SDSS collaboration who made our observing runs and the SDSS 
project a success. 
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Fig. 1. — The W A2796 distribution. Left: The distribution of rest equivalent widths, W^ 2796 , 
for Mg II systems found in the survey. Right: The redshift-path covered by the survey, 
AZ (W X2796 ) = f* max Ef spec ^W 2796 , z) dz, as a function of W* 2796 . 
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Fig. 2. — The redshift-path corrected Wq 2796 distribution from the MMT data. The mean 
absorber redshift is (z) = 0.589. The solid line is a maximum likelihood fit of the form 
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to data having W^ 2796 > 0.3k, with W* = 0.511 ± 0.047 and 



N* = 1.071 ± 0.119. The red-dashed line is the low-z result from the SDSS EDR survey 
with (z) = 0.655. The points represent data from CRCV99 with (z) = 0.9. The MMT result 
confirms the upturn in ON/OWq 2796 below 0.3A that was first identified in NTR05. 
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Fig. 3. — Evolution of W*. The points without symbols are from the EDR survey and the 
circles are the MMT result for Wq 2796 > 0.5A. The curve is the power law fit described in 
NTR05. 
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Fig. 4. — The absorption redshift distribution. Left: The distribution of redshifts for Mg II 
absorption systems found in the MMT survey. Right: The total number of sightlines with 
sufficient signal to noise ratio to detect lines with W A2796 > W^ in as a function of redshift, 
for W^ in = 1.0 A, 0.6A, and 0.3A. 
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Fig. 5. — dN/dz as a function of look-back time, for a (fl\, Qm, h) = (0.7, 0.3, 0.7) cosmology. 
The horizontal bars represent bin size and the vertical bars the la uncertainties. Points 
without symbols are the EDR results, while the circles represent the MMT and x symbols 
the CRCV99 results. The dashed lines are the no-evolution curves normalized to the binned 
EDR/CRCV99 data, but excluding the MMT data. 



